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The cubic Complex Ginzburg-Landau Equation (CGLE) has a one parameter family of traveling 
localized source solutions. These so called 'Nozaki-Bekki holes' are (dynamically) stable in some 
parameter range, but always structually unstable: A perturbation of the equation in general leads 
to a (positive or negative) monotonic acceleration or an oscillation of the holes. This confirms that 
the cubic CGLE has an inner symmetry. As a consequence small perturbations change some of the 
qualitative dynamics of the cubic CGLE and enhance or suppress spatio-temporal intermittency in 
some parameter range. An analytic stability analysis of holes in the cubic CGLE and a semianalytical 
treatment of the acceleration instability in the perturbed equation is performed by using matching 
. and perturbation methods. Furthermore we treat the asymptotic hole-shock interaction. The 

' results, which can be obtained fully analytically in the nonlinear Schroedinger limit, are also used 

for the quantitative description of modulated solutions made up of periodic arrangements of traveling 
holes and shocks. 
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I. INTRODUCTION 

> 

Spatially extended oscillatory media in the weakly nonlinear regime can be described by a complex Ginzburg Landau 
equation (CGLE) which may be derived from an amplitude and slow- variation expansion Q. In Id it reads 



d t A =[! + (! + ib)d 2 x - (1 + ic)\A\ 2 + d\A\ 4 ] A (1) 



o 

. 

where b and c are real constants . Length, time and the complex order parameter A have been scaled in the usual way 
for the cubic CGLE. In cq.(@) we have included a quintic term with a small complex prefactor d = d! + id" (\d\ <C 1). 
We will refer to eq.(Q) with d ^ as the "perturbed (cubic) CGLE" and treat the d term only perturbatively to lowest 
. order 0(d). 

As already pointed out in reference Q arbitrarily small \d\ may change the situation drastically, and eq.(^) without 
I- fifth order term (d = 0) is not sufficient to describe the qualitative dynamics of the system in the parameter regime 
considered here, i.e. where Nozaki-Bekki hole solutions (see below) influence the dynamics. From the derivation of the 
• >"H , CGLE one expects in a real physical system d cx e where e is the distance from threshold. There are other corrections 
to the cubic CGLE which should have the same order of magnitude but their " perturbative effect" is expected to be 
rather similar and we have chosen the d term in eq. ([!]) to represent " higher-order perturbations" . While in reference 
we presented a phenomenological theory for phenomena observed numerically, in this work we will determine the 
important quantities of the theory with the help of analytical or semi-analytical methods. 
Equation (|l|) has the family of traveling-wave solutions 

A = R q e l< - qx - u ^ t] with R q = ^l-q 2 + d'{l-q 2 ) 2 , < q 2 < 1 (2) 
and the dispersion relation 

w(q) = cR\ + bq 2 - d"R\ . (3) 
(As mentioned above all formulas are valid only at order d .) 



A. Localized Solutions 



Other important solutions are localized objects A v moving with constant velocity v connecting asymptotic plane- 
wave states. In the co-moving frame (C = x — vt) the plane-wave dispersion relation reads fl(q,v) = io{q) — vq. 
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Dealing with such constantly moving solutions it is convenient to write the CGLE (|lj) in a co-moving frame rotating 
with constant frequency f2 



d t A(C,t) = F v , n [A] 

:= (1 + in + (1 + ib)d 2 c + vd c - (1 + ic)\A\ 2 + d\A\ 4 ) A (4) 

Clearly, in order to ensure phase conservation, it is necessary that in this co-moving frame the plane waves on both 
sides rotate with the same frequency, i.e. Lo{qi) — vqi = O for i = 1,2. From this and eqs.(^) and (||) one can derive 
the relation: 

2(b-c)k-2(cd'-d")(Rl+R 2 q Jk (5) 
with ft = + q 2 ) (6) 

Note that in the cubic CGLE this yields v — 2(6 — c)k which is the mean of the group velocity d q u>(q) of the 
asymptotic plane waves. One differentiates between sources and sinks (often called shocks) depending on whether the 
group velocity (and thus causality) points inward or outward in the co-moving frame. Thus sources should determine 
the behavior of the asymptotic states while shocks are passive and play only a minor role. So in the following we will 
focus mainly on the investigation of sources. 

For a standing (v = 0) (antisymmetric) source solution eq.(|J) is fulfilled trivially due to symmetry (qi — — (ft)- For 
v =/= 0, however, eq.(|^) is not trivial at all and, as pointed out by van Saarloos and Hohenberg the problem to 
find such source solutions of eq.(^) is in general overdetermined. Using counting arguments they showed that from a 
systematic point of view only standing sources should exist. Nevertheless the cubic CGLE (eq.(l) with d = 0) has an 
analytic one-parameter family of sources with different velocities v, the "hole solutions" of Nozaki and Bekki |Q (NB 
holes). They can be written in the form 

A^ B := Z v exp(iXv ~ iSU) 

= [BdQ(p v (n() + Av} x exp[iLp v (n() + ikv( — i£lt] (7) 
where (p v {nC) = k^ 1 lncosh(/c£) . 

Symbols with a "hat" denote constants depending only on b and c, e.g. k = 1/(2(6— c)) . The frequency fl and k 2 are 
linear functions of v 2 . With exception of A all coefficients are real. The asymptotic plane waves for ( — * ±oo have 
wavenumbers 

qi/2 = k ± K with K = k/ k , k = vk . (8) 

k 2 becomes zero at a maximal velocity ±v maxi and here the hole solution merges with a plane wave with wavenumber 
<7i = q2 = v max /(2(b—cj). The relations are completely derived in Appendix The resulting algebraic equations (8 
equations for 8 parameters) yield the one parameter family. The derivation in Appendix |A| demonstrates the necessary 
dependence in a transparent way. 

The NB holes are connected with the 3-parameter family of dark soliton solutions that exist in the limit b = 
±oo, c = ±oo , where eq.([lj) reduces to the defocussing nonlinear Schroedinger equation (see next paper), and with 
the one-parameter family of static (all velocities collapse to zero) saddle-point solutions of the nonlinear diffusion 
equation obtained in the limit b = c = Then one has a potential and the continuous family can be associated 
with (global) gauge invariance through Noether's theorem. Clearly then perturbations of the equation that preserve 
gauge invariance and the potential property, like a quintic term with a real prefactor, do not destroy the family. 

An important result of this article is that in general the family of hole solutions is destroyed by the higher-order 
perturbation (d ^ 0) leaving only the standing hole. Thus moving NB holes are structurally unstable. 



aj{q 2 ) - c4gi) 

qi - qi 



B. Stability and Destruction of Hole Solutions 



The stability of the hole solutions (for d = 0) has first been investigated by Sakaguchi in direct simulations of 
the CGLE. The perturbational equation of the CGLE (Q) with A = A v + W describing the stability problem 

d t W = C v W with C V W := -^[A V ]W + -j^[A v ]W* (9) 
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(where now A v — A^ B and SI = 0(t?)) describing the stability problem was then studied by Chate and Manneville 
numerically for v = 0, v = 0.2 ||] and by Sasa and Iwamoto semianalytically for v = [0]. As a result hole solutions 
were found to be stable in a narrow region of the b-c plane which is shown in Fig.l for the standing hole (v = 0) . 
From below the region is bounded by the border of (absolute) stability of the emitted plane waves with wavenumber 
q(b, c) (see eq.(^)) corresponding to the continuous spectrum of C v . From the other sides the stable range is bounded 
by the instability of the core with respect to localized eigenmodes corresponding to a discrete spectrum of C v . The 
core instability turns out to be connected with a stationary bifurcation where the destabilizing mode at threshold 
passes through the neutral mode W = Qtr (see eq.(|l4|)) of C v which can be derived from translational invariance of 
the CGLE ^||. The cores of moving holes were found to be more stable than those of the standing ones 

Now we turn to the situation with small but finite d. Our simulations show that in the stable range moving holes 
are then in general either accelerated and eventually destroyed or slowed down and stopped to the standing hole 
solution depending on the phase of d. In particular for real d — d! one has 



d t v 



v 



> 0^d' > (10) 

< < 



One finds that the Nozaki-Bekki hole relations connecting the core velocity and the emitted wavenumbers (see eq. (0) 
and App.|A|) are (almost) satisfied at each instant during the acceleration process. The acceleration thus occurs 
approximately along the NB hole family and it can be described by taking v = v(t) as a slowly varying variable while 
other degrees of freedom follow adiabatically. 

From these numerical observations we conclude that for arbitrarily fixed values of b,c and varying d the standing hole 
solution undergoes a symmetry-breaking stationary bifurcation at d = 0. The (real) growth rate A^ of the unstable 
mode is proportional to the acceleration 

Ad = lim v/v . (11) 

We will refer to this instability as the " acceleration instability" while the term " core instability" will be used for the 
bifurcation where b and/or c is changed while d is kept constant (including the case d= 0). 

Of special interest is the case where the core-stability line is crossed with d (which actually corresponds to the 
physically relevant situation). Then the two modes which cause the acceleration instability and the (stationary) core 
instability (with d = 0) are coupled which leads to a Hopf bifurcation, (nur bei d stabilisierend) As we showed in j|] 
the normal form for this bifurcation - valid for small d, v, u — is 

u = (\ — sv 2 )u + div . . 

v = fXU + d2V . 

Here u and v are the amplitudes of the core-instability and acceleration-instability modes respectively, s and [i are of 
order d° while d\$ must be of order d 1 since in the absence of a perturbation holes with nonzero velocity exist. For 
d = the parameter A can be identified with the growth rate of the core-instability mode. (In the general case A has 
a correction proportional to d causing a shift of the threshold.) The nonlinear term in cq.([l2"|) takes care of the fact 
that moving holes are more stable than standing ones and at the same time saturates the instability (s > 0). 

Far away from the core-instability threshold where A is strongly negative u can be eliminated adiabatically from 
eqs.(p^) which yields 

v = (d 2 -d 1 ^)v (13) 



The term in brackets can be identified as A^ from eq. (|1 1|) . 

After introducing our general tools and concepts in Secjl| we will present a detailed perturbation analysis in Sec. Ill 



Since the core instability (for the case d = 0) as well as the acceleration instability are of the stationary type (real 



growth rate) we will first investigate the stationary bifurcation scenario more closely in Sec. [II A. Our approach 



allows to treat both instabilities within the framework of the same formalism. Then we will derive a fully analytic 



expression for the core instability line in the unperturbed equation (d = 0) (Sec. Ill B), present a semi-analytic method 
for treating the acceleration instability in the case d ^ (Sec. |lII C| ) and investigate the interaction between holes and 



shocks in Sec. Ill 3 . A comparison of our analysis with simulations of the CGLE is presented in Sec. IV. Finally we 



will use our results for a quantitative description of moving periodically modulated solutions of the perturbed CGLE 
(Sec.0). 
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II. TOOLS AND CONCEPTS 



A. The Linear Operator 

The linear operator C v introduced in eq.(^) has two bounded neutral modes 

<P rot = iA v <P tr = d c A v (14) 

related to gauge and translational invariance of the CGLE. 

For many applications it is convenient to introduce a transformation W — W exp(— i$ v ) where is some real field 
(often the phase of A v ). Equation (||) then goes over into 

d t W = L V W := exp(-ii?„)£ t , (exp(i$ v )W) . (15) 

For the NB solutions (d — 0) we choose $ v — \v obtaining 

L% B W = ((1 + ifi) + (1 + ib)(d c + id cXv f + v(d c + id cXv ) ~ 2(1 + ic)\A* B \ 2 ) W-(l + ic)\A^ B \ 2 W* 
= (l+ib)d cc W + (2Ki{l+ib)tanh(K0+2ki(l+ic))d i W 

+ ((l+ic)(-K 2 B 2 tanh 2 «) + 2fcXtanh«) - 2v 2 A 2 ) + (1+ifi + (l + ib)iKn)/ cosh 2 ( K C) + ivk - k 2 (l + ib)) W 



- ((l+ic)(vA + KBt&nh(K()) 2 ) W* . (16) 

Stationary instabilities are described by the eigenvalue problem related to eq.([l5|) 

L V W = XW (17) 

with real A. The bounded neutral modes from eq.(|l4|) are now 

$ rat =iZ v = i(vA + A'Btanh(KC)) (18) 
$ tr = d ( Z v +iZ v d c Xv = KBk/ cosh 2 «) + (k + Ktanh(n()) $ rot . (19) 



The other two fundamental modes satisfying L^ B ^> = (irrespective of boundary conditions) can be expressed in 
terms of generalized hypergeometric functions. 

Asymptotically 3> 1) one has plane waves, the operator L v becomes space independent and the four funda- 
mental modes behave exponentially ~ e Pi< ^ . The exponents pi (i — 1...4) can be obtained from the characteristic 
polynomial. Two of the exponents can be extracted from the bounded neutral modes $ ro t and <&tr They are 

Pl = and p 2 — t2k for ( — > ±oo (20) 

The other - in general complex - exponents p$ and p4 are given in Appendix |b|. Their real part is always posi- 
tive/negative for C < showing that the other two fundamental modes are exponentially growing in space. 



B. Condition for the Existence of a Family 

A necessary condition for a standing hole solution Aq := A v= q of the CGLE to be embedded in a continuous 
one-parameter family of moving hole solutions A v can be obtained by taking the derivative of the stationary CGLE 
-Fufi[AJ = with respect to v at v — : 



0«7ctm 



toa) [A ] = dM =#tr (21) 

dA v I 

f am dv \v=0 



with <£ ■ - 



(Here the fact has been used that for symmetry reasons fl can only depend on v 2 .) Equation ( |2l| ) shows that the 
existence of a hole family implies that the translation mode <S>tr = d^Ao has an inverse image ^/ am under Cq := £ v =o 
which for large |£| grows linearly in space (^f am ~ iC-A v ) corresponding to an asymptotic wavenumber change (c.f. 
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eq.(|2^)). In the cubic case (d — 0) eq.(pl|) with this boundary condition is solved by inserting the Nozaki-Bekki 
solution A v = A^ B . Then instead of eq. (f2l|) one may write (c.f. eqs.QqMph) 



Ln B & 



jam 



with - <f> 



fam 



A + ikBKC,t-<mh(KC,) 



(22) 



In general (d ^ 0) all solutions of eq. (^l|) are found to diverge exponentially for £ — > ±oo (see below) which is consistent 



b-c 



l+bc 



<C 1 , which 



with the fact that the hole family is destroyed by the higher-order perturbation. In the limit ^ 
includes the relaxative case b = c = as well as the conservative and fully integrable nonlinear Schroedinger limit 
(b, c — > oo), we show this analytically in ||. 



C. Motion within the Family — Asymptotic Matching 



Above we stated that the acceleration of holes in the perturbed CGLE can be (approximately) described as a 
motion within the family of Nozaki-Bekki hole solutions by taking the family parameter v — v(t) as slow variable. 
This would, however, lead to a change of the asymptotic wavenumbers (one has k — kv) which is not possible in an 
infinite system. (As pointed out before the family mode <£/ m = d v A v {(,,t) , eqs.( pl|p2[ ), describing the difference 
between two "neighboring" holes diverges linearly for large \(\ .) The difficulty is resolved by limiting the "motion 
within the family" to a finite region of size ~ v^ 1 around the hole core (inner region). A global solution for A has 
then to be constructed by asymptotic matching. For calculating the acceleration, however, it will not be necessary to 
perform the matching explicitly. 

To illustrate this we consider the eigenvalue equation ( p"7| ) with the condition 

\W\ < bounded for £ -► ±oo (23) 

describing the acceleration of the hole in the limit of small velocities (see next section). In the outer region 

|«C| > 1 (outer region) (24) 

W behaves exponentially with exponents pi = Pi(X) (c.f. Sec. pl A| , App.|B|). For values < A <C 1 there are two 
decaying (pi/2) and two growing (P3/4) exponents so that boundary condition ( p3| ) is indeed equivalent to 

'W should not grow in space exponentially (with the exponents P3.4)' . (25) 

It is crucial to note that because of the overlap of regions the realization of this boundary condition can already be 
controlled in the inner region without performing any matching. 

Since there is one weakly decaying exponent p\ ~ A ~ v boundary condition (|25] ) in practice leads to an algebraic 
growth of W (at the outer limit of the inner region). This can be understood as a truncated expansion of the 
exponentially decaying outer solution (see eq.([Bg|)) in the overlap region W ~ (i + 0(A)) (f + pi( + 0(X 2 )) . The 
portion ~ ipiC corresponds to the linear divergence of $j am . From this one finds the outer limit of the inner region 

|piC| "C 1 (inner region) . (26) 

Boundary condition (p5| ) is not restricted to the case of small velocities. Applied to the inner region it remains also 
valid in more realistic situations with slowly varying wavenumbers governed by phase equations far away from the 
hole core. We note that besides perturbations of the equation also core instability or modified boundary conditions 
(coming i.e. from the presence of a shock) in general lead to a hole acceleration which can be treated by an analogous 
matching approach. 



III. PERTURBATION ANALYSIS 



A. Formal Analysis of Stationary Instabilities 



In this section we will investigate stationary instabilities of the standing hole with respect to localized modes. The 
formalism will be applicable to the acceleration instability in the perturbed CGLE and to the core instability in the 
cubic CGLE (d = 0). 
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Let u = (b, c, d) denote a point in parameter space of the CGLE . We write u = u + Au where u belongs to the 
stability threshold. In particular u may lie on the core- instability line (b cr it,c cr it,d = 0), then Au = (A6, Ac, 0) is 
of interest (core instability), or u can be a point (b, c, 0) where the standing hole is stable and then Au = (0,0, d) 
(acceleration instability). We want to construct the unstable localized mode W(x; u = u + Au) near the threshold. 
It satisfies 

L Q (u)W(u) = A(u)W(u) (27) 

where Lq has been introduced in eqs.(|s|JT5|). 

In the following we want to expand W(u + Au) in powers of Au. The calculations will be limited to the inner 
region (eq.^6j) where one may use boundary condition ( |25| ) instead of the usual condition ( |23| ) for localized modes. 
We will use the symbol Au to denote the distance from threshold along a fixed path which crosses the neutral curve. 
Consistently the derivative signifies the derivative along this path. At the threshold, A(u) = 0, we have 

W(u) = $ tr (u) . (28) 

To see this we note that, since ^(u) belongs to the null space of £o( u ), it has to coincide with one of the two neutral 
modes given in eq. (|l4|) . For reasons of symmetry one expects to have <I>t r , which destroys the zero at the hole core. 
Inserting the expansion 



W{u + Au) = $ tr (u + Au) + AuFTi(u) + Au 2 W 2 (u) + 0(Au 3 ) 

<9A(u) 



into eq.(27) one finds at first order in Au 

L (u)Wi(u) 



Equation ( pfj| ) is proportional to eq. (f2l|) and thus has the solution 

<9A(u) 



Wi 



du 



4> 



farn 



u) 



At second order in Au eq.(^) then yields 

id 2 A(u), ... 



L (u) W 2 (u 



d\(u) f dXju) 

du \ du 



$/™(u) 



a$ tr (u) 9L (u) 



(29) 



(30) 



(31) 



(32) 



where the term 
with 



1 o 2 a(u) 

2 du* 



$t r (u) was brought to the lhs with the help of eqs.([ 



du du 

fj) and (HI). If a function $/ am (u + Au) 



L (u + Au)$ /am (u + Au) = $ tr (u + Au) . (33) 

exists at u + Au the term in square brackets can be absorbed by modifying W2 ■ To see this one may expand W in 
the following way 

dA(u) 



W(u + Au) = $ tr (u + Au) 
which yields at second order in Au : 



du 



-$ /am (u + Au) + Au z W 2 (u) + O(Au^) 



L (u) ( W 2 (u) - - 



1 d 2 A(u) 



du 2 



jam 



(u) = 



d\{u) 

du 



$/am(u) 



(34) 



(35) 



Equation (p3[), which is the condition for the existence of a hole family at u + Au, is satisfied near the core-instability 
line if d is kept zero. Thus the solvability of eq.([55|) gives a criterion for the threshold of the core instability (d = 0) 
which will be evaluated fully analytically in the next section. Near the threshold of the acceleration instability however 
eq.(|3^) can only be fulfilled if 

<9A(u) _ d 



du 



dd 



(36) 



d=0 



is adjusted properly which can be used to determine the acceleration of the holes. Thus the quantitative results for 
the acceleration of the holes given below represent a proof for the destruction of the hole family by the higher-order 
perturbation. 

From eq.(^2|) one also sees that v/v diverges as one approaches the core-instability line in the b-c plane. Near this 
line however the analysis of this chapter looses its validity. Here the two modes leading to core and acceleration 
instability interact and the scenario becomes more complex (see eq. (fL2|) ) . 



G 



B. Core Instability 



In this section we calculate the core instability line in the cubic CGLE. We do this by giving a solvability condition 
for eq. ( j35|) which we write in the form 

L(u)W(u) =$/ ow (u) . (37) 

Here $/ am is the family mode from eq.(p2|) and L(u) = Lj^ B (6, c) is the linear operator around a standing hole (see 
eq.dlfl)). Equation (|37|) reads explicitly 



(l+ib)d a W + (2iK(l+ib)tanh( K 0) d c W + (1 + iQ + (l+ib)iK K ) Tl — -W - 

cosh [kQ 

- (l+ic)(l-/sT 2 )tanh 2 (KC) (W + W*) = -K - 2ni + i(l-K 2 )Ctavh(KC) (38) 

Solvability here means that we have to adjust the parameters b — 6(c) in a way that the (symmetric) solution W of 
this equation satisfies boundary condition ( p5| ) . In a first step we simplify eq. (|3^) by making the following ansatz for 
W ("variation of constants" ) 

W = f rot® rot + ftr^tr 

where $ rot = i tanh(/«£), $ tr = 1/ cos1i 2 (k£) + (i/k) tanh 2 («;£) (39) 



are the two neutral modes of Lq (c.f. eqs.( [L8|]19| )). f ro t and f tr are two real functions to be determined. Inserting this 
ansatz into eq.(38) one gets after elimination of f ro t the following real second-order differential equation for / := d^ftr 

Lf = If am 

where L := <9 CC - 6Ktanh(/0<% + ((6K 2 + 12k 2 ) tanh 2 (/<) - 4k 2 ) 
I fam ■- 6^ 2 Ccosh 2 ( K C) + 3K 2 (bk - 2)C + (b + bK 2 - 4K 2 k)3K/(l - K 2 ) taah«) cosh 2 (/<) (40) 



From the asymptotic behavior {\(\ — > oo) of cq. (J40[) one easily sees that the only boundary condition for the (anti- 
symmetric) function / = c^(cosh 2 (nC,)Re(W)) which is compatible with condition ( |25| ) for W is given by 

/ - Ce 2 " ICI for C -> ±oo . (41) 

This excludes the behavior ~ e ( 2K +P3/4)KI Q f the ^ w0 fundamental modes of L. Equations (40 4l]) can be solved with 
the ansatz 

oo oo 

f~ J2 c m Ccosh" 2m (KC) + ^mtanh(KC)cosh- 2m «) (42) 

m= — 1 m— — 1 

As shown in Appendix |^ for certain values of (6, c) with b = 6(c) one can adjust the coefficients c m , d m in a way 
that the series converges uniformly towards a smooth function solving the boundary value problem eqs.( [40| , f4ll) . The 
condition for convergence defines the core instability line. Explicitly this line is given by 

^S(k 2 ) + 1 + K 2 ^ = (43) 



_,„ 2 , 3 ^ 2n+l -Pj 2m 2 -m + 3 k 2 

where S(k 2 ) := — - + > — — - '— (44) 

V ' 4 + 3 k 2 ^ 2n 2 + 2n + 3 k 2 l±2m 2 + m + 3 k 2 K ' 

' 71—2 ' m — 2 ' 

_ 3 ^ r(|+a)r(|-a) g 2n+ 1 I> + f + a)I> + | - a) 



4 + 3/k 2 T(| + a)T(| - a) ^ 2n 2 + 2n + 3/k 2 T(n + § + a)I> + | - a) 



with a := - J\ - 24/k 2 



4 

In the limit of large c one can neglect the first term in eq.(f43|) which contains the infinite sum S(k 2 ) and eq.(|4£ 
reduces to 
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either b 4 = ^c 2 + 0(c 3/2 ) for c -» oo (b > 0) (45) 
9 

or & = -l/v / 2 + 0(l/c) forc->oo (46) 

Eqs.(^) describe both branches of the core instability line in the limit c — > oo. The core instability line (|43|) together 
with the curve describing the instability of the plane waves at the wings of the hole give the complete stability diagram 
of the standing Nozaki-Bekki hole solution. The curves are plotted in Fig.l. They are consistent with previous results 
jjj |§ U and generalize them. Finally we remind that our derivation made use of the fact that the core instability 
occurs via a static bifurcation, which was first found by B and H. 



C. Acceleration Instability 

In the following we treat the acceleration of a moving hole caused by a perturbation of the cubic CGLE. We use 



adiabatic elimination which in the limit v — > can be put into the framework introduced in section III A . We assume 
that in the perturbed CGLE the hole acceleration occurs near the hole family and can be described by taking v = v(t) 
as the slow variable while the other degrees of freedom follow adiabatically. Formally this is done by writing 

A=(Z v(t) (C) + W(C)) e X p(^ Xv{ t)(0-ij\^v(r)+^)dry, ( := x - J\(t)<It (47) 

where Z v ,\v, &v describe the slowly time dependent NB hole. W, AT2 are the small changes of the solution caused by 
the quintic perturbation. At lowest order they are time independent (adiabatic elimination). As before this ansatz is 
only valid in the inner region Q. Inserting @ into the CGLE (g) in the co-moving frame rotating with frequency 
£l v (t\ + Afi one obtains 



vd v A = F v(t)Qv(t)+m [A] (48) 



The only time dependence which cannot be transformed away by this choice of the coordinate system is given by 
a movement within the family vd v A. W, v, AT2 are of order d. Using the fact that A^ B solves the cubic CGLE 
F^^[A^ B ] = at each time t and neglecting all higher order terms in d one arrives at the following ordinary linear 
differential equation for the perturbation W 

L N V B W = I v (Ml,v) 

where I v (Ml,v) := u$/ om - i£flZ v - d\Z v fZ v (49) 



Here one has to adjust the parameters AT2, v in the inhomogeneity in such a way that the solution W satisfies (|2£ 
which determines the acceleration it. 

Equations ( |32| ) and ([l9]) can both serve to calculate the acceleration sufficiently far away from the core instability 
line (where it is possible to replace eq.(|l2|) by eq.(|l3|)). While cq.(f49|) is applicable for holes moving with arbitrary 
velocities v (but small v) eq.(|3^) is only valid for small v. The connection between both equations can be seen by 
expanding ( pj| ) in terms of v. Writing 

L» B = L» B W+vL% B W+0(v 2 ) 
W = W {1) +vW {2) +0(v 2 ) 

I v = + vl {2) + 0(v 2 ) (50) 

and remembering that L^ B ^ = Lq B is the linear operator from eq.(|38|) describing perturbations around the standing 
hole of the cubic CGLE, eq.(|49|) becomes at order v° 

L$ B W^=I^ 

= -imz - dz\ (51) 

Here 1^ and thus W w are antisymmetric. has to satisfy the boundary condition ( p5| ) which fixes the parameter 

and Afi describe the changes of a standing hole resulting from a d-perturbation in the cubic CGLE. These 



changes have already been included in the approach of subsection III A. and Afi can be found analytically (see 



cq.(D13)). At order v 1 eq.(ff9[) then reads 
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= -L*[ B MWW + -$/am " «AT!i - dZ*(3A + 2A*) (52) 



v 



This equation corresponds to eq.(|32|). The r.h.s., and thereby the perturbation W^ 2 \ are symmetric and one has 
to adjust the acceleration - such that fulfills the boundary condition (^5|). (Note that in deriving eq.(^) we 

have expanded the equation of motion (0) in d and afterwards in v whereas eq. (^) can be found from eq. ( Eg) by an 
interchange of the expansions.) For small d,v one has 



- =Re(g*d) + 0{v 2 ) g:= 9l +ig 2 (ft real) (53) 



In Section IV A eq.([i"9|) will be used to determine g numerically for arbitrary b, c while in App.Jo] eq.(|5^) will be used 



to calculate g analytically in the NLS-limit (b, c— > oo) 



D. Hole — Shock Interaction 



The acceleration of a hole as calculated in the last subsection is influenced in the presence of a shock. Here we treat 
this acceleration in the limit of large hole-shock separation. For that purpose we divide space into overlapping hole 
and shock regions. The important point is that when using the ansatz (^) in the hole region the boundary condition 
( p5|) is changed by the shock. In the region of overlap one now has (besides decaying ~ e Pl ? 2 '') growing stationary 
perturbations ~ e Pi/^ forming the shock out of the plane-wave state. The distance from the shock determines the 
magnitude of the growing perturbation (prefactor of W) and thus the hole acceleration. 

The hole-shock interaction depends crucially on the fact whether the growing exponents p 3 / 4 are real or complex 
conjugated. Accordingly one has monotonic ('monotonic case') or oscillatory ('oscillatory case') interaction (see below 
eqs.( |57|j5^ )). The boundary between the two interaction types in the b, c-parameter plane is given by the condition 
Pz = pi which depends on the hole velocity v since p 3 / 4 = £3/4(1;) (c.f.App.|B]). In particular for a standing hole the 
boundary line becomes 

k 2 — 6 = (54) 



which is included in Fig.l (re is given in eq.(A8).). 

In the overlap range (re -1 <C C ^ L) between the hole position £ — and the shock position £ = L (defined by 
d^arg(A)\^ = L — 0) one has approximately plane waves plus linear perturbations W of the form (c.f.AppjB|) 

We-** = r 3 e P3C+ ^ 3 + r 4 e P4C+#4 (monotonic case) (55) 
or We~ ilp = ze P3<: + 7y(p 3 )z*e p ^ (oscillatory case). (56) 

(Here <f> :— arg{vA + KB) is a fixed phase factor.) The perturbations are made up of the two growing fundamental 
modes of L v . The two real prefactors r 3 / 4 = r 3 / 4 (L) or the complex constant z = z(L) describing the contribution 
of these modes have to be determined from the asymptotic shape of the shock solution (see below). They depend 
exponentially on the distance L between hole and shock with reversed exponents ~ e~ Pi ' iL . Thus the L dependence 
of the acceleration induced by the perturbation W is given by 

v = ,93e~ P3L + g 4 e~ P4i (monotonic case) (57) 
or v — g 3 e^ Ret - P3)L sin(Im(p 3 )L + g^) (oscillatory case) . (58) 



The two unknown real constants 33,34 have to be determined from solving eq. (|49|) with boundary condition 



instead of (25). In the monotonic case the smaller exponent p 4 (p 3 > p 4 > 0) dominates asymptotically and for 
practical purposes the coefficient g 4 is (usually) sufficient to describe the acceleration. 

Now we turn to the calculation of the constants r 3 / 4 or z from the asymptotics of the shock solution. Besides the 
long-wavelength approximation, where after a Hopf Cole transformation the shock region is described by a linear 
phase equation (see M), there exists no analytic expression for the shock solution. Therefore in general the constants 



r 3 ,r 4 ,z have to be calculated numerically by solving a boundary value problem (cf.Sec.1V A). In addition a crude 
analytic estimate can be obtained by the following 'linear approximation': we extend the hole region up to the center 
of the shock which means that the shock is a superposition of a plane wave and the linear perturbations ~ e P3 / 4< ». 
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Assuming d^A\Q = L = at the position £ — L of the shock (which corresponds to a standing shock) one gets by using 
again @ 

9 C ((^ ( *) + M/ ) elx " <t) )l c =L =° 
(d f W + % (vA + KB + W^ ^ =0 (59) 

Inserting W from eq.(|55|) or (^) one obtains an approximate expression for r^/^L) or z{L) respectively. 

So far we have only treated the interaction of one hole and one shock in the cubic CGLE. In general boundary 
conditions ( |55| , p6| ) apply on both sides of the hole separately with corresponding shock distances L( r /iy In most cases, 
however, one of the hole-shock interactions dominates and the constants r 3 / 4 or z can be set zero on the other side 
(corresponding to the interaction with a shock infinitely far away). In any case the full information on the boundary 
condition for W can be gathered into one real four vector b containing the boundary constants r 3 or z for both 
sides. Since in the perturbed CGLE the hole acceleration depends linearly on this boundary vector b and on the d- 
perturbation (via the inhomogeneity of eqs.( |49|j5^ )) one may simply add the hole accelerations caused by neighboring 
shocks (eqs.( |57[|5"8| )) and d-perturbations of the CGLE (eq.(|5^)). This is exploited in the next section. 

IV. NUMERICAL RESULTS AND COMPARISON WITH SIMULATIONS 

A. Numerical Methods 

We now show how one can determine numerically the hole acceleration induced by a d-perturbation or an interaction 
with a shock. Mathematically this means solving the differential equation ([l9|) with the boundary conditions (|25| ) or 
(|55|,^6|) described by the boundary vector b. 

First let us assume the vector b were known. One has to adjust v, A£l in eq.([49|) in a way that the solution W 
is asymptotically described by b. This can be achieved by the following 'shooting' method. One integrates eq.(|49|) 
(numerically) starting at £ = 0. The correct initial values W(0) and d(W{Q) are unknown but one may assume without 
loss of generality Re(W(Q)) = and Im(d^W(0)) = 0, which can always be achieved by a suitable addition of the 
neutral modes <& ro t and $tr to the perturbation W. One is left with two real parameters Im(W(0)) and Re(drW(0)) 
describing the initial condition and two further real parameters v, Afi occuring in the inhomogeneity of eq.(^9|). Now 
one integrates eq.(|4^) from £ = once to (sufficiently large) positive £ and once to negative £ with all four parameters 
set zero. From the asymptotic growth of W at the wings one can read off a real four vector bo describing this growth. 
Doing similar integrations with the following values of the parameters one obtains the boundary vectors hi, (i = 1...4). 

Im(W(0)) = 1 Re{d c W(0)) = v = AT2 = =S> b l 

Im(W(0)) = Re(d c W(0)) = 1 i = AT2 = b 2 , , 

Im(W(0)) = Re(d c W(0)) =0i = im = 0^b 3 ^ ' 

Im(W(0)) = Re(d ( W(0)) = w = Ml ^ 1 ^ fe 4 

Because of the linearity of the boundary value problem one can now determine the set of parameters corresponding 
to the desired boundary vector b. Solving the linear equation 

b Q + 7m(W(0))(6i - bo) + Re(d c W(0))(b 2 ~ &o) + v(b 3 - b ) + Afi(6 4 -b ) = b (61) 

yields the acceleration v and the frequency shift AT2 , which the hole receives under the influence of neighboring shocks 
and/or a d-perturbation. One needs two of these runs for the calculation of 31,52 (eq.([53|)) and two further runs for 
the calculation of 53,54 (eqs.( |57|j5^ )). 

Now we turn to the determination of the constants r 3 ,r^z in eqs.( |55| , p6| ). Since the shock solution connects two 
plane-wave states it is described by the ODE (|j) F^ = ^ [A] = with the shock velocity v s and frequency fl s fixed by 
the (given) wavenumbers q\, q 2 at the wings. We solved this boundary value problem using a NAG routine. In some 
intermediate region between the system boundaries and the center of the shock the amplitude is just given by the 
incoming plane waves and linear perturbations ~ e P3/iC growing towards the shock center. In this region we extract 
the desired constants r 3 , r 4 , z. 

Our simulations of the CGLE (Q) were performed with a pseudospectral code based on FFT using a predictor 
corrector scheme in time. We had to use very high precision in the simulations, because discretization errors, like 
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the d-perturbation, in general destroy the inner symmetry of the cubic CGLE and therefore have similar effects on 
the hole solution, i.e. they also lead to an acceleration. Typically we used discretizations of Ax ~ 0.05 — 0.15 and 
At ~ 0.001 — 0.02 depending on the parameters 6, c, d of the CGLE. 



B. Results 



Next we compare the results from the semianalytical calculations in sections III C III D with those of direct simu- 
lations of the CGLE. 

Figure 2 shows the acceleration of a hole center caused by a quintic perturbation of the CGLE. The acceleration 
formula fl5S| ) with the constant g obtained from the numerical matching (using eqs.(^9 25)) is seen to be in good 
agreement with the simulations. This also confirms that the acceleration instability is indeed of stationary type. 
When scanning the (b, c)-range of stable hole solutions one always find s dece leration for (real) d < 0, consistent with 
the simulations. From the fully analytical treatment in App.|j5| (c.f. eq.(D14)) it is seen that this result remains valid 
in the nonlinear Schroedinger limit. (Then one has numerical problems with the matching because here the ratio 
\Ps/Pa\ of the two growing exponents P3,Pi tends to infinity which prevents a correct numerical extraction of the 
boundary vectors hi (c.f. section IV A ).) 

As already stated in section III D the (b, c, cQ-region of stable hole solutions is divided into parts with monotonic 
and with oscillatory interaction between holes and shocks. Figures 3a, 3b show the acceleration of a standing hole 
under the influence of a standing shock solution for both interaction types in the case d = 0. Here in the direct 
simulations the shock solutions were simply modeled by taking d^A = at the system boundaries. As a result one 
has almost perfect quantitative agreement of simulations and numerical matching if in the latter the shock is treated 
nonlinearly, i.e. when the parameters (or z) in the boundary conditions ( |55| , [56|) a re extracted from a numerical 

solution of the shock boundary value problem. With the 'linear approximation' (|59|) one obtains only qualitative 
agreement. When scanning the monotonic interaction range with the matching routine we found that the asymptotic 
hole shock interaction is always attractive. Near the nonlinear Schroedinger limit one has monotonic interaction (see 
Fig.l). Here one can calculate analytically the hole-shock interaction (see |^]) and it also turns out to be always 
attractive. 



V. ARRANGEMENTS OF HOLES AND SHOCKS 

In this section we use the results from the previous sections for the description of stable states made up of a 
periodic arrangements of holes and shocks. Such states are frequently observed in simulations with periodic boundary 
conditions (see e.g. Figure 4 shows the modulus \A\ = |A(x)| of a typical solution found in a simulation. 

As shown in Fig. 5 one finds uniform as well as (almost) harmonic and strongly anharmonic oscillating hole velocities. 
(Slightly) beyond the core instability line the direction of the velocity is changed in the oscillations (see Fig.5d). The 
solutions are seen to be very sensitive to d-perturbations of the cubic CGLE. 

The uniformly moving solutions can be well understood from the results of the last sections. First one has to note 
that they are not expected to be stable (and we indeed could not observe them in si mulat ions) in the monotonic range 
[] since here the asymptotic hole-shock interaction is always attractive (c.f. section fVE). Therefore we may restrict 
ourselves to the case of oscillatory interaction. Away from the core instability line uniformly moving periodically 
modulated solutions can then be identified as fixed points of a first order differential equation for the hole velocity v 
as the only slow variable. The other degrees of freedom are thereby assumed to follow adiabatically. 

One equation for the hole velocity v is obtained by considering the acceleration caused by interaction with a 
neighboring shock at distance L together with that resulting from a d-perturbation. From eqs.(^3|,[58|) one obtains in 
the limit of small velocities v 



vRe(g*d) + g 3 e-^ L sm(pUv)L + g 4 ) , 



9 = 9i+ W2 



P3 



P3 



ip'l 



(62) 



Here we have taken into account only the interaction of the hole with one of the neighboring shocks, implicitly 
assuming that the period P of the solution satisfies P 3> L in addition to kL ^> 1 , which is necessary for the 



The boundary between monotonic and oscillatory interaction depends on the hole velocity (see HID) 
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asymptotic analysis. Equation ( |62[ ) is correct up to first order in v if the two matching parameters #3,4 are calculated 
including corrections linear in v and L in the matching procedure (i.e. g^iv, L) = 930 + vg^ v + Lg 3 j^). In the following 
we eliminate the two unknowns L, L. 

The shock velocity v s = v + L is determined by the incoming (plane) waves via a generalization of eq. (|^) to slowly 
varying solutions leading to 



v s — v + 0(v) 



(63) 



Thus L — 0(v) can be neglected in the matching parameters 53,94 occurring in eq.(|62|) (since the factor e~ P3 ^ L is 
small) . 

From periodicity one has the condition that the phase along one period P of the solution is an integral multiple of 
2ir. Using the fact that the full phase of a NB-hole solution is given by 



4>(0 = Xv(C) + arg(Z v (0) 

w if ICI -\n2 + k( + arg(vA±KB) for ± k( > 1 



one can derive the phase condition 



2nn = (j>(P) -0(0) w 2KL + P(k - K) + arg 



vA + KB 
vA-KB 



+ <p s {v) (n e N) 



(64) 



(65) 



Here 4> s is the phase change due to the presence of the shock. For symmetry reasons cj) s — <p s (91,92) with qi = qi(v) 
is an antisymmetric function of the hole velocity v and it becomes linear for small velocities. Solving (|65| ) for L and 
linearizing in v one gets 



L*Lr. 



g^v 



with: L n p 



(P + (2n-l)ir/K), g 5 :-- 




P 



u=0 



2A" 
KB 



Inserting eqs.(|63|) and ( |66[) into eq. (|62|) one arrives at a differential equation for v of the desired form 

V = f n p(v) 



(66) 



(67) 



describing solutions with n, P as parameters in addition to 6, c, d. In the derivation we have assumed stationary (or 
at most slowly evolving) solutions as well as small velocities. Furthermore we have neglected cZ-corrections in the 
plane wave quantities (k, K,p3^, ...). All parameters (gi,-.,95) can be extracted from the matching procedure and 
the numerical nonlinear shock solution. 

The implications of eq.(|6^) are very simple. The fixed points (v = 0) describe the uniformly traveling states. Their 
velocity Vf p = Vf p (n,P) and hole-shock separation Lf p = Lf p (n,P) are predicted by eq.(|67|) as functions of n,P. 
Their stability (within the family parameter v) is determined by the sign of the linear growth rate A/ p := df n p/dv\ Vfp . 
In Figs. 6, 7 these theoretical predictions for Vf p (n, P) and L f p {n, P) are compared with data from numerical simulations 
for the CGLE parameters b = 0.5, c = 2.3, d = ±0.0025. The full lines describe stable solutions (A/ p < 0) while 
dashed lines correspond to unstable ones (A/ p > 0). The theoretical results are seen to be well confirmed by the 
numerical simulations especially in the limit of small velocities. The remaining discrepancies can be explained by the 
neglect of the d-corrections in the plane wave quantities (fc, K,p^^, ...) and the linearization in v. Note that for fixed 
period P several stable solutions can coexist which may (but don't have to) differ by the phase 27m contained in one 
period (see Fig. 6). 

The essential influence of the e?-perturbation on the solutions is especially apparent for very slowly solutions (c.f. 
eq.^f 



v fp « -(1/Re(g*d)) 930 e~ p '^ L ^ sin(p 3 '(0)L„ P + 940) 



(68) 



occuring in the limit of small interaction e~ Ps " p <C \d\ <C 1. Since one then has A/ p w Re(g*d) these solutions are 
only stable for decelerating c£-perturbation. As a consequence for accelerating d (slowly) uniformly moving solutions 
exist stably only for certain periods (see Fig. 7) whereas in the opposite case such solutions occur for arbitrary values 
of P (see Fig. 6). This fact is well confirmed in our simulations. 

Now we turn to the solutions with oscillating hole velocities. In our simulations we found such solutions coexisting 
(stably) with uniformly moving solutions but they occurred especially for parameters where eq.(p7|) does not have 
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stable fixed points (see Figs. 5b, 5c). The solutions are characterized by the occurrence of large velocities (v ~ 1) 
and interaction strengths (e Pa ( v ' L ~ 1) during the oscillations. Presumably they can (in a first approximation) be 
identified as stable limit cycles of a two dimensional dynamical system with the hole velocity v and the hole-shock 
distance L as active variables. While one has L = 0(v) (as in eq.(|63|)) v, and thus L, cannot be neglected now. Thus 
from the phase conservation condition one obtains an evolution equation for L which has to be combined with the 
equation for the hole acceleration v yielding the two dimensional dynamical system. In this description cqs.( p6p7| ) 
correspond to the adiabatic elimination of L. 

The mechanism for the oscillations beyond the core stability boundary (see Fig.5d) is of different nature. Here 
the core instability mode introduces a new degree of freedom. While for accelerating d this leads to a destruction of 
the solutions, in the decelerating case the interplay of core and acceleration instability leads to a Hopf bifurcation 
of the fixed point solutions of eq.(|67[) similar to the situation described for isolated holes by eq. (p"2|) . This yields the 
oscillating solutions of Fig.5d. 



VI. CONCLUDING REMARKS 



The investigation presented in this paper is of general importance since the CGLE is generally derived from an 
amplitude expansion near the threshold of a bifurcation and higher-order perturbations appear naturally. Of course 
there should be other terms of the same magnitude as the quintic term included in eq. (Jl|) . However, since the 
effects described above are closely related to phase conservation one should expect the other terms to have the same 
qualitative consequences and it should be possible to describe a physical system with a single "effective" higher-order 
term. We expect the consequences of our results to be restricted to (quasi-) 1-dimensional systems in the vicinity of 
the parameter range where hole solutions are stable in the (unperturbed) cubic CGLE. Thus stable, uniformly moving 
holes should in general not be observable. (Note that in the case where the background state carries traveling waves 
rather than homogeneous oscillations, the CGLE ([!]) has an additional group velocity term. Motion then has to be 
defined with respect to a frame moving with that group velocity.) 

The detailed investigation of spatio-temporal chaos presented by Shraiman et al. for the unperturbed CGLE is 
presumably robust with respect to perturbations because the sensitive parameter range is not considered there. A later 
investigation jll|] considers the sensitive range, where the determination of the range of existence of spatio-temporal 
chaos becomes questionable in the unperturbed CGLE. We found that stable hole solutions suppress spatio-temporal 
chaos and as a consequence for a stabilizing perturbation the (upper) boundary of spatio-temporal chaos is simply 
given by the stability boundary HS (see Fig.l) of the NB hole solutions, whereas for destabilizing perturbation 
spatio-temporal chaos is also observed further up . 

Comparing holes in Id with spirals in 2d one can state that spirals behave in a rough sense (apart from topological 
stability) similar to holes in the perturbed CGLE with stabilizing (decelerating) perturbation. Thus one has a range 
(which is in 2d very large) where standing spirals are core stable. For large values of |&| a core instability occurs, 
which is of Hopf type [[L3| . In contrast to the Id case the bifurcation is subcritical in 2d and leads to spatio-temporal 
chaos. We also note that the theory of interaction between localized solutions is similar in Id and 2d [H . 

Transient hole-type solutions were observed experimentally by Lega et al. |T^ 1 in the (secondary) oscillatory in- 
stability in Rayleigh-Benard convection in an annular geometry. Here one is in a parameter range where holes are 
unstable in the cubic CGLE, so that small perturbations are irrelevant. Long time stable stationary holes ('Id spirals') 
were observed in a quasi-ld chemical reaction system (CIMA reaction) undergoing a Hopf bifurcation by Perraud et 
al. |lq| . The experiments were performed in the vicinity of the cross-over (codimension-2 point) from the (spatially 
homogeneous) Hopf bifurcation to the spatially periodic, stationary Turing instability. Simulations of a reaction- 
diffusion system (Brusselator) with appropriately chosen parameters exhibited the hole solutions (and in addition 
more complicated localized solutions with the Turing pattern appearing in the core region). It seems likely that for 
the hole solutions the vicinity of the codimension-2 point is not important, and that the parameters correspond to 
the stable hole range in the CGLE description with stabilizing perturbations. 

Finally we mention recent experiments by Leweke and Provencal JlTj where the CGLE is used to describe results 
of open-flow experiments on the transitions in the wake of a bluff body in an annular geometry. Here the sensitive 
parameter range is reached and in the observed amplitude turbulent states holes should play an important role. 

We wish to thank S.Sasa for sending us data of the core instability line in addition to those in ||. One of us (LA.) 
wishes to thank the Alexander- von-Humboldt Stiftung for financial support and the University of Bayreuth for its hos- 
pitality. Support by the Deutsche Forschungsgemeinschaft (Kr-690/4, Schwerpunkt 'Strukturbildung in dissipativen 
kontinuierlichen Systemen: Experiment und Theorie im quantitativen Vergleich') is gratefully acknowledged. 
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APPENDIX A: THE NOZAKI-BEKKI HOLE SOLUTIONS 



Inserting the ansatz (0) into the CGLE one obtains a polynomial of third order in dQip v {nC : ) = K tanh(K^) (K := 
k/k). Requiring that terms at each order cancel separately one is confronted with 4 complex equation for 8 real 
parameters so that a priori one should only expect a discrete solution. The 4 complex equations, however, are not 
independent, which leads to the one parameter family of hole solutions. In the following we derive the explicit form 
of the Nozaki-Bekki hole solutions in a way which shows this dependence from the structure of the CGLE and the 
solutions. (Results are listed in eqs. (|A| - |Al|)). Introducing (@) into the CGLE one obtains 



(Lu n [v, k, fi] - (1 + ic)\Z v \ 2 ) Z v e 1 ^ = , Z v = B0 c¥ >„«) + Av 



(Al) 



with the linear operator Li in [v, k, Q] := 1 — k 2 + i(Q + v 2 k — bk 2 ) + V8q + (1 + ib)d 2 
and V :— v(l + 2i(l + ib)k = vV). Using the fact that Lu n commutes with the operator (—iBd^ +vA) one easily finds 



L Hn Z v e 1 ^ = L Hn (-iBd c + vA)e^ = 



AZ V - Bid c A 



v» 



(A2) 



with the "eigenvalue" A = A(d^ip v ) :— exp(—iip v )Lu n exp(iif v ) = z(d^ip v ) 2 —iVd^^p v + l+i(n+v 2 k)+i(l+ib)(kK 2 — k 2 ) 
and z := (1 + ib)(ik — 1). Introducing this into eq.(Al) one sees that d^A/Z v has to be a polynomial of second order 
in d^ip v = if tannic) and one finds 



A = zB~ 2 Z 2 



const. 



from which follows 

iV = 2zAB- x . 
Using eqs.(A2) and (A3) equation ( |A1[ ) reduces to 

A - (1 + ic)\Z v \ 2 = 2izB~ 1 d c Z v 
which is of second order in d^(p v . Dividing this equation by z the real part has the form: 

.1 



B~ z Re{Z 2 ) ~ Re(- 



-)\Z V \ = const. 



(A3) 



(A4) 



(A5) 



(A6) 



(where the fact that d^Z v — Bd^<p v is real was used). Noting that Re(Z 2 ) = (d^ip v + vA') 2 — A" 2 = \Z V \ 2 + const one 
finds that the equations obtained from the nonconstant contributions - those proportional to (d^ify) 1 and (c^,,) 2 - 
are equivalent. In this way one can convince oneself that the remaining equations are indeed not independent. 

To get the explicit solutions from eq.(A5) one may proceed as follows. At order (d^(p v ) 2 one obtains a complex 
equation 



(1 + ic)B 2 = (1 + 2ik)z = (1 + 2m)(l + ib)(ik - 1) 



which determines B and k: 



and 



4(6 -c 



^- (3(1 + be) + V9(l + &c) 2 +8(6-c) 2 ) 



B 



3k(1 + 6 2 ) 

b-c 



(A7) 



(A8) 



(A9) 



From order d^tpl contributions of eq.( A5) one gets another complex equation which together with eq.( A4) yields a set 
of 3 independent real equations for k,A' and A". Solving them one obtains 



A' = -B- 1 A" = and k = [2(6 - c) 

B 



(A10) 
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where the last formula is identical with the phase-conservation condition (|5|). From constant terms (~ (d^ip v )°) one 
can derive the dispersion relation (in the moving frame) 

n = c-vk + (b-c)(k 2 + K 2 ) (All) 

and an elliptic relation which relates the two not yet specified variables v and n (or K = k/k) to each other: 

v 2 (k 2 + \A\ 2 ) + K 2 (1 + B 2 ) = 1 . (A12) 

Thus one finds indeed a one parameter family of moving hole solutions which can be labeled e.g. by their velocity. 

-l 



From eq.( A12) one finds that k becomes zero at v = v n 



APPENDIX B: ASYMPTOTIC ANALYSIS OF NOZAKI— BEKKI HOLE SOLUTIONS 



In the limit ( — > ±00 the NB hole solutions (7|) become plane waves with wavenumber q 1 / 2 '■= k ± K (sec eq.(@)) 
and frequency ^(^1/2) — c + {b— c )Qi/ 2 ( see ec f(|3))- The perturbational equation for these plane waves in a coordinate 
system moving with the velocity v of the hole core reads 

d t W = (1+ib) d cc W + (v + 2iq 1/2 (l+ib)) d c W - (l+ic)(l- q\ /2 ) (W + W*) (Bl) 

This equation coincides with the asymptotic (\k(\ S> 1) perturbational equation (|l5| ) with L v from eq. (|l6|) (after 



transforming away an (irrelevant) fixed phase factor W -> We- iar 9( vA±KB )). Equation mm can be solved by an 
exponential ansatz 



with p, A, z, complex 



(B2) 



Inserting this into eq.(Bl) yields two complex linear equations for the constants z, z* whose solvability condition leads 
to characteristic equation 



p\l + b 2 ) + p 3 2v + p 2 (±q 2 1/2 + (v- 2bq 1/2 ) 2 - 2(1 + bc)(l - q\ /2 ) - 2A) 

+ p ((4(6 - c)q 1/2 - 2v)(l - qf /2 ) + (Abq 1/2 - 2v)\j + (a 2 + 2A(1 - <? 2 /2 )) 







(B3) 



The four roots pi = Pi(v, A) describe the possible asymptotic spatial behavior of perturbations W with growth rate 
A in the coordinate system co-moving with the velocity v of the hole core. In particular the exponents pi = piiv) 
describing stationary perturbations are obtained for A = 0. For £ — > +00 one finds 



Pi =0, 

P3/4 = K 



v 



-2k 



a 2 02,0/1 2\ 1 + bc , ^bqiv + 2nv - v 2 

YTb 2 ± V^- 3k2 + 2 ( 1 -^tt^ + 



1 + 1 



{i + b 2 f 



For v — one has q\ = —q 2 and the expression for p 3 / 4 simplifies 



P3/4 = K± Jk 2 - 6q 



(B4) 



(B5) 



The exponents p% for C — > —00 can be found from eqs.(B4,B5|) by the replacement q\ —> q 2 and reversal of the 
overall signs. In the following we will for simplicity only treat the case £ — > +00. The exponents p\ and p 2 describe 
the asymptotics of the two neutral modes eqs.([l8],[l9]). The exponents P3/4 being real (monotonic case) or complex 
conjugated (oscillatory case) always satisfy 



Re(p 3/4 ) > 



(B6) 



which shows that two fundamental modes of L^ B are exponentially growing in space. For each pi the constants 
z(pi), z„(j>i) can be found from the eigenvector equation. One finds 



15 



= rjz 

For real pi (which covers the exponents pi,P2 and in the monotonic case also P3-,Pi) one has = 1 and the 

asymptotic behavior of the stationary perturbation associated with one of the exponents pi is given by 

w „ e P.C+#. where 0. := -arg(rj(pi)) for C -> oo { Pl real) (B8) 

In the case p 3 = p% complex W behaves asymptotically like 

W = ze PsC + r){p 3 )(ze P3C )* for ( -> oo (p 3 complex) (B9) 

The complex quantity z gives the two growing fundamental modes in the perturbation W . We note that depending 
on the velocity v of the hole the behavior of the growing perturbations ~ e P3/4.C m ay be quite different on both wings; 
in particular it may be oscillatory on one and monotonic on the other wing. 



Now we turn to the situation of a small nonzero growth rate (|A| -C 1). Then the exponents from eq.(B4) are slightly 
changed. In particular the exponent p\ = pi(v,X) describing the outer asymptotics of localized modes bifurcating 
through the translation mode $t r (e.g. the core instability mode) becomes near their bifurcation 

Pi~ ^-At vT fOT l A l«! ( B1 °) 

v - 2 1 gi (6 - c)| 

Performing a similar analysis of asymptotic plane wave states in the perturbed CGLE (0 ^ \d\ <C 1) most of the 
foregoing expressions receive d-corrections. However, since translational (and rotational) invariance are preserved by 
the perturbed equation, the stationary (A = 0) exp onent p\ — remains unchanged and also eq. (|B10| ) for small A 



remains valid at lowest order in d. Therefore eq.( BlO ) describes (besides the core instability mode in the cubic CGLE) 



also the acceleration mode in the perturbed CGLE. 

APPENDIX C: CORE INSTABILITY 

In this section we present the solution of the boundary value problem eqs.(^(],|d]) defining the core instability line 
utilizing the ansatz (fl2|). From the relations 

L (Ccosh~ 2 "(KC)) = coslT 2rl «) (XC - M£Ccosh~ 2 «) - i?„tanh«)) (CI) 

L (tanh(KC)coslT 2 ™(KC)) = cosIi~ 2 "(kC) (7V n tanh«) - tanh(/«C) cos1i~ 2 (kC)) (C2) 

with N n := 4n 2 K 2 + Ylnn 2 + 8k 2 + 6K 2 , R n :=4Kn + 6K (C3) 

Ml := 4n 2 K 2 + Unn 2 + 12k 2 + 6K 2 , := 4nV + 18/jk 2 + 20k 2 + 6K 2 , (C4) 

follows, that when inserting (^) into eq.(|40|), one can in a first step recursively determine the coefficients c m and 
then in the next step the coefficients d m from the c m . In summary one has 

c-i = 1, co = ^rp, (C5) 

= ^, (CT) 

*. - + - ±- (ft fe) { + + t k l 3) } CO 



where 



Ai := 2k 2 + 3bnK, A 2 := j^(3b + 3bK 2 - UkK 2 + 2k) (C9) 



1G 



From eq.(C6) follows that for large m the coefficients c m behave like 



Ar_ 'iff Ml 1_ 1+jk 
N m 11 N n m2 11 1 + 1 



m 3 ' 2 for rn — * oo 

I 4- -i 

n=0 n=0 

Therefore the first sum in (143) always converges uniformly towards a smooth function. Because of 



(CIO) 



1 

"Mr, 



m—1 

n 

71=0 



N„ 



d 



1 Tr 1 ! 



TT 2 " 



n=0 



- 1 / 2 for 



(Cll) 



the coefficients d m in eq.(C8) may behave like ~ m~ 3 ' or ~ m ' depending on whether the curly bracket in 
eq.(C8), which behaves like const + 0(l/m) , gives zero in the limit m — > oo or not. In the latter case the second 
sum in the ansatz eq. (^2|) converges (pointwise) towards a function with a finite step at C, = whereas in the first case 
the limiting function is smooth everywhere. Therefore one has to adjust (6, c) with b — b(c) in a way that the curly 
bracket in eq.(C8) vanishes for m — ► oo. This determines the core instability line. Resubstitution of N m , M^, ... 
into the curly bracket of eq.(C8) yields the final expression (|43|). 



APPENDIX D: ACCELERATION INSTABILITY IN THE NONLINEAR SCHROEDINGER LIMIT 



In this section we present a method to treat the acceleration instability in the NLS-limit (6, c — > oo) by making use of 
the formalism introduced in subsections II C , III A . IIIC . We calculate the acceleration by projecting the (symmetric) 
inhomogeneity of eq.([52|) onto the symmetric neutral mode Q s of the adjoint operator ^ := Lq B . The scalar 
product vanishes only for suitably adjusted parameter - in the inhomogeneity and this ensures the solvability of 
eq S .(p||). 

Using the scalar product 



the adjoint problem is 



(U,W) 



Re(U*W)dx 



L^ s = 



with 



14, 



for C — > ±oo 



(Dl) 



(D2) 



(l-ife)<9 cc $ s + (2iir(l-i6)tanh«)) d^ s + (1 - + (l-ifyiKn) 



1 



cosh 



4., 



- (l-ic)(l-if 2 )tanh 2 (KC) $ s - (l + ic)(l-^ 2 )tanh 2 «) $* 



In analogy to Appendix M an investigation of the asymptotics of the fundamental modes of shows that for ^> 1 
they behave like ~ e Pi '^T For C — > oo the exponents pt are 



Pi = 0, 



p 2 = +2k, 



P3/4 



-K ± \/h 



6K 2 



(D3) 



These are just the exponents of L with reversed sign (c.f. eqs.(B4,B5)). The two decaying exponents P3/4 describe 
the desired symmetric mode <fr s in the outer region (kC\ S> 1) of the asymptotic matching procedure. In the NLS 
limit one finds with the help of the relations of App.jA 



P3/4 



-re ± k 1 



4(6 



36 2 c 2 



0(6- 4 ) 



(D4) 



Both exponents are real ('monotonic range') and ps — > for b,c — > 00. The contribution to $ s proportional to e P3 ^ 
has a prefactor Z3 with 



2 In this section we drop the indices '0' and 'NB' to simplify the notation. 
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Im(z 3 ) 
Re(z 3 ) 



+ 0(b' 3 ) 



(D5) 



as can be calculated in analogy to (B8). The inner region is defined by [psC\ — \ K C ^%^ I ^ 1 which ensures a finite 



overlap in the NLS limit. Here we expand the differential equation ( p2[ ) in terms of 1/6 

L = b(L<® +b- 1 L^ +0(b- 2 )) 



^ =b[L<M +b- 1 L^ +0(b- 2 ) 
$ s = $i 0) +6- 1 $i 1) +O(6- 2 ) 



At order b 1 eq.flDSj) reads 



(2tanh 2 (0 ~ l)*i 0) ~ -tanh 2 (/<)$i 0) * = 



(D6) 



(D7) 



from which one gets $ s ° ' = i/ cosh 2 (K£). At order b° eq.(D2) yields 
50 . L (o)t $ (i) = _ L (i)t$(o) 

= i (tanh 2 «)(8 - 28k 2 ) + (-4 + 8k 2 )) /(3 cos1i 2 (kC)) 

Using the full fundamental system of L^ ^ (irrespective of boundary conditions) 

ij cos1i 2 (k£) 

tanh(«^) l — k£ tanh(K£) 

one can construct the bounded solution 

^ = ^^tanh 2 K) + l 
3k z 



(D8) 



^- cosh(KC) sinh«) + tanh(KC) + ^ J 

4k 8k 8 cosh (kC) 



(D9) 



(D10) 



Combining $1°'' and and taking the limit C — ► 00 in the terms 0(b 1 ) one finds the outer behavior of the inner 
solution 



$ s = ;/cosh 2 (KC) + + Q(b- 2 ) 



(Dll) 



This expression has to be matched to the outer exponential solution. By using eqs.(D^-D5,Dll) one finally gets 



cosh (kQ 



2 -4k 2 
36k 2 



wICI h_i l +0 (b- 3 )) (l + (9(6- 1 )) 



(D12) 



Now we turn to the inhomogeneity of eq.(^). For arbitrary 6, c the terms AO, occuring here can be found 

from the following ansatz describing the standing hole of the perturbed cubic CGLE in the inner region 



tanh(/xia;) + /13 tanh 3 ^!^)) exp (i^ In cosh.(fiix) + i/15 tanh (/iix) — ifigt) with fa real 



(D13) 



Taking the NLS limit of the inhomogeneity in eq.(^2|) and projecting it onto the neutral mode <!> s then yields the 
acceleration 



v ( 36 c 



8 6 26 



v \2(b-c) 3 



1 + 0(b-')) +----_ + 0(6- J ) ) = -Re(d) 
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15 



(D14) 



As stated in Sec. Ill A the acceleration v/v diverges at the boundary of the core instability (eq.(f45|)). 
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Figure Captions 



Figure 1: Stability diagram for the ID CGLE. Waves emitted by the standing hole become convectively unstable 
below the dashed curve EH (Eckhaus instability) and absolutely unstable below the solid curve HS. Outside the 
region bounded by the solid curve CS the core of the standing hole becomes unstable (Hopf bifurcation for (real) 
d < 0). The CS line is given by eq.(f43|). Stable standing holes can be found above HS and below CS (upper 
branch) for d < 0. The dashed dotted line MO gives the boundary between monotonic (above) and oscillatory 
(below) interaction for standing holes (see eq.(p4)). 



Figure 2: Acceleration instability in the perturbed (cubic) CGLE. Comparison of the reduced acceleration dtv/v 
from theory (full line) and simulations (squares) for b = 0.5, c = 2.0, |d| = 0.002 and varying phase arg(d). 
The theoretical curve was found from adjusting (numerically) the parameters v and AQ in eq.([49|) leading to 
v/v = i?e(0.6572 e 66901 d), see eq.Q. 

Figure 3: Acceleration v of a standing hole in the presence of a shock at distance L in the cubic CGLE (d = 0). 
Comparison of simulations (squares) and theory. The parameters r 3 , r<4, z in the boundary condition eqs.( p4|]55| ) 
were obtained from the (numerical) nonlinear shock solution (full lines) or from the analytic approximation 
eq.(|5rj) (dashed lines), (a) monotonic range: b = 1.0, c = 3.5; full line: v = 3.395 e~ PiL , dashed line: v — 



4.255 e~ PiL with Pi = k - y/n 2 ~ 6q 2 « 0.4877 (see eq.@). (b) oscillatory range: b = 0.5, c = 2.3; full 



line: v = 4.361 e'^^ 3 ^ sm(Im(p 3 )L - 1.460); dashed line: v = 8.823 e ~ Re( - p ^ L sin(Im(p 3 )L - 0.874); with 
p 3 = K + ^k 2 - 6q 2 w 0.8896 + 0.5940i (see eq.@). 

Figure 4: Snapshot of the modulus \A\ — \A(x)\ of a stable uniformly moving hole-shock pair in a simulation for 
b = 0.5, c = 2.3, d = +0.0025. The solution is space periodic with period P = 48.4. 

Figure 5: Simulations showing the velocity of the hole center v = v(t) of interacting hole-shock pairs (periodic 
boundary conditions). In Figs.(a),(b),(c) the CGLE parameters were b = 0.5, c = 2.3, d = +0.0025 (i.e. far 
away from the core instability line), (a) relaxation into a constantly moving solution for period P — 48.4. (b) 
selected final state with (almost) harmonic oscillating velocity for period P — 37.0. (c) selected final state with 
anharmonic oscillating velocity for period P = 40.0. (d) CGLE parameters b — 0.21, c = 1.3, d = —0.005 (near 
the core instability line); state with oscillating velocity (including change of the direction) for period P = 50. 

Figure 6: Interacting hole-shock pairs. Velocity v — v(n,P) and bound state distance L = L(n,P) are plotted as 
function of the period P for different phases 2mr contained in one period for the CGLE parameters b = 0.5, c = 
2.3, d = —0.0025. Full lines show stable (A/ p < 0) and dashed lines unstable (Xf p > 0) uniformly moving 
solutions which are predicted (fixed points) by eq.(|67|). They are compared with (stable) uniformly moving 
solutions observed in simulations (X). For larger velocities there are more stationary and oscillatory solutions, 
which have not been included since the analysis is not applicable in that range. 

Figure 7: Like Fig. 6, but with CGLE parameters b = 0.5, c = 2.3, d = +0.0025. 
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